Electron Removal Self Energy and Its Application to Ca2Cu02Cl2 
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We propose using the self energy denned for the electron removal Green's function. Starting 
from the electron removal Green's function, we obtained expressions for the removal self energy 
Y, ER (\c,u) that are applicable for non-quasiparticle photoemission spectral functions from a single 
band system. Our method does not assume momentum independence and produces the self energy 
in the full k-cj space. The method is applied to the angle resolved photoemission from Ca2Cu02Cl2 
and the result is found to be compatible with the self energy value from the peak width of sharp 
features. The self energy is found to be only weakly k-dependent. In addition, the luiE shows a 
maximum at around 1 eV where the high energy kink is located. 

PACS numbers: 74.72.-h, 74.25.Jb, 78.20.Bh, 79.60.-i 



I. INTRODUCTION 

Photoemission lineshape contains information on the 
many-body interactions in the solid under study. One 
notable example is the kink structure in the angle re- 
solved photoemission(ARPES) spectral function from a 
metal at the bosonic mode energy due to electron-bosonic 
mode couplingi^. In this regard, the kink structures 
at near 70 meV found in the ARPES spectral func- 
tions from cuprate high temperature superconductors 
(HTSCs) have been a hot topic during the last several 
years as it may provide a clue on what mediates the pair- 
ing in the HTSC^&1£. It is still under debate whether 
the kink is caused by phonons or magnetic excitations. In 
addition to the kink near 70meV, an anomalous high en- 
ergy kink structure from insulating cuprate Ca2Cu02Cl2 
(CCOC) was recently reported^. Subsequently, it was 
found that the high energy kink also exists in doped cases, 
making it omnipresent at all dopings 10,11 . Existence of 
the high energy kink reveals that there may be another 
energy scale from the many-body interactions in HTSCs. 
It is thus important to find the physics in this high energy 
scale. 

One can extract from the data the information on the 
many-body interactions through the self energy analy- 
sis. Conventionally, the self energy analysis is done by 
extracting the dispersion and peak width from the mo- 
mentum distribution curves (MDCs) 1 . However, there 
are underlying assumptions for the analysis to be valid. 
First of all, ARPES peaks should be sharp enough to 
be quasiparticle-like. Secondly, the bare band dispersion 
should be quite linear within the energy range of interest 
so that one can "guess" the bare band dispersion. Un- 
fortunately, none of these are valid assumptions in the 
analysis of ARPES spectral function from CCOC: The 
peaks are too broad to be considered as quasi-particle 
peak a 12 i 13 and there is an inherent confusion regarding 
what should be considered as the original bare band as 
CCOC is a Mott insulator. Similar limitations also apply 



to the recently found high energy kinks in HTSC a 10 ^ 1 . 
However, the lack of quasiparticles should not be con- 
fused with a lack of physics. For example, in Luttinger 
liquids, where Fermi liquid quasiparticles do not exist, 
ARPES has provided detailed information on the excita- 
tion spectrumJ^. 

It is therefore necessary to develop a method to extract 
the self energy from non-quasi particle ARPES spectral 
functions. In fact, there was an earlier attempt to obtain 
the self energy through a Hilbert transformation with the 
assumption of particle-hole symmetr}^. This method 
however is not applicable to systems without particle- 
hole symmetry (for example, CCOC—). Here we propose 
to use the self energy defined for the electron removal 
Green's function. Even though different, it bares similar 
characters of the usual self energy and its imaginary part 
is found to appropriately represent the width of the spec- 
tral feature. Moreover, we developed a method to extract 
the self energy from a general ARPES spectrum from a 
single band system. It can thus be used to for spectral 
functions with broad spectral features. We believe our 
analysis method provides a new way to examine the un- 
derlying physics from photoemission spectral functions. 



II. METHODOLOGY 

The full Green's function for a single band system at 
K can be written as 
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where Sk and E Fn "(k,c<;) are the bare band and self en- 
ergy, respectively. The chemical potential \i is set to 
zero for convenience. Note that due to the sum rule 
E n {K n M0)| 2 + |(n|4|0)| 2 }=l, the asymptotic behav- 
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ior (uo — > oo ) of G Ful1 is l/o;. This makes the self energy 
non- divergent when uo — > ±00. In addition, the self en- 
ergy is analytic in the positive complex plane, and as a 
consequence the real and imaginary parts of the self en- 
ergy are Hilbert transforms of each other—. The Green's 
function defined in this way ("full Green's function") is 
thought to represent electron addition or removal pro- 
cesses and the self energy holds the information on the 
interaction. 

The full Green's function in the above equation how- 
ever includes not only electron removal but also electron 
addition process for which we do not have information. 
Our goal is to obtain the self energy from a general spec- 
tral function with as few assumptions as possible. There- 
fore, we attempt to look at the self energy for the electron 
removal Green's function which should best describe the 
ARPES process and is defined as 
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where n k = ^Z n \(n\c k \0)\ 2 = J A ER (k, uo)duo and 
A ER (k,uo) is the electron removal spectral function (the 
superscript ER stands for "electron removal"). Another 
difference is that e k may not be the bare band dispersion 
as will be discussed later. The asymptotic behavior of G 
is now n k /uo and with the definition given in equation 
(2), the real and imaginary parts of the electron removal 
self energy (ReE FjR and ImE^, respectively) are again 
related through Hilbert transforms. 

To obtain the electron removal self energy E FjR , we 
use the fact that the measured ARPES spectrum is pro- 
portional to the electron removal spectral function, that 
is, 

Iarpes(^)=I = |M(k)| 2 A^(k,^) 

= -M^ ImG - (k ^) (4) 
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where Iarpes and M(k) are the measured ARPES spec- 
trum and matrix element, respectively. Here, the uo de- 
pendence of the matrix element M(k) is ignored as usual. 
Note that the Fermi function is missing as we deal with 
the electron removal Green's and spectral functions. The 
real and imaginary parts of the electron removal Green's 
function are, as is the case for G Ful1 , Hilbert transforms 
of each other, KeG ER =H(lmG ER ) where H() represents 
the Hilbert transformation. Then the combination of 
equations (3) and (4) givesi^, 
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where e k and C(k) are the bare band and -n k \M(k)\ 2 /ir, 
respectively. 

The only remaining task to get E is to determine C(k) 
and e k . We can obtain C(k) by integrating both sides of 
equation (4). 
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To define e kl we enforce the standard relationship 
ReE £;i? =H(ImE £;i? ). For this relationship to hold, 
Y> ER (k,uo) — > as uo — > ±00. By expanding the last term 
in equation (3) in terms of 1 /uo for uo — > ±00, we get 
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Since Yi ER (k,uo) — > as uo — > ±00, we set e k = e k and 
thus we have 
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f A ER (k,uo)duo ~ J I ARPE s(k,uo)duo 



(10) 



That is, e k is the center of gravity of the spectral function 
but may not be the true bare band. 

We rationalize this definition of e k from the knowl- 
edge that e k = the bare band value, indepen- 
dent of the strength of interactions in the limit n k — > 
Alternatively, e k could have been chosen differently 
such as the true bare band value of the original non- 
interacting problem but it will only shift ReE^ deter- 
mined here by a constant. In that case, the relation- 
ship between ReE FjR and ImE^ would be modified to 
ReZ ER =R(ImZ ER )+e k -e k . We note then that ReE^ 
goes to a finite value of e k -e k as uo — > ±00 (similarly 
ReE Fn " asymptotically goes to the Hartree-Fock value). 
Using this expression we can see that the parameter e k 
drops out and that the imaginary part of the self energy 
remains uniquely determined independent of our choice 
of e k . 



III. RESULTS 

In applying the newly developed method to the 
ARPES data, one has to consider other aspects of the 
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FIG. 1: (Color online) (a) Published ARPES data from 
CCOC along the (0,0) to (tt,tt) cut in Ref 8. (b) Same data 
after background is removed. Data is presented for a smaller 
momentum window where peaks are reasonably strong. 



experimental data such as experimental resolution and 
thermal broadening. In that respect, CCOC can be an 
ideal system to test the method: The experimental res- 
olution and thermal broadening are much smaller than 
the width of the characteristic features of CCOC and 
thus can be ignored. This provides an extra motivation 
for this study. For the experimental data, we use the 
published CCOC data along the (0,0) to (7r,7r) cut shown 
in Fig. 1(a) 9 . 

Before we apply the procedure to extract T, ER , we also 
have to subtract the background intensity from the data. 
This is because the background intensity is from a pro- 
cess that can not be described by the Green's function 
in equation (2). In our case, subtracting the background 
so that A ER (k,uj) — > as u — > ±oo is a necessary con- 
dition for ReE^ and ImE ER to be Hilbert transforms 
of each other. The background intensity is also casually 
subtracted in other methods, either naturally by fitting 
the MDCs^ or enforcing h— . There are a number of po- 
tential ways to subtract the background, among which we 
may consider the following two. The first one is subtract- 
ing the Shirley backgroun d 19 ^ . The other is subtracting 
the intensity of the unoccupied k states (approximately 
beyond k = 0.7A in Fig. 1(a)), that is, the MDC back- 
ground assuming that it is from the momentum indepen- 
dent background intensity. We show the result extracted 
by using the latter method in Fig. 1(b). However, the 
final results are similar independent of the background 
subtraction method. 

In panels (a) and (b) of Fig. 2, we show ReE^ and 
l m ^ ER obtained by applying the procedure discussed 
above. Note that while ReE^ can be negative, -ImE^ 
is always positive. The most notable aspect of the data is 
that overall k dependence for both ReE FjR and ImE^ is 
weak. This is very striking considering the fact that the 
dispersion of the spectral features is very strong. The 
momentum independence of the self energy which has 
been casually assumed ^ 22 1 23 1 24 ^ is indeed seen in the 
data. The momentum independence could further be 
tested for our systems. Yet, we still see some momen- 
tum dependence of the self energy, especially at the low 
energy side. Therefore, for the accurate measure of the 
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FIG. 2: (Color online) (a) ReE Si? extracted from the data in 
Fig. 1(a). The energy scale is coded in the color scales on the 
right, (b) -lmE ER . The line is the dispersion extracted from 
fitting the MDCs. (c) -lm^ ER along the MDC dispersion in 
panel (b). The inset shows the EDC from the k point where 
the spectrum is the peakiest. The full width at half maximum 
of the peak is about 300 meV. (d) Comparison of the Sk (solid 
line), and the MDC dispersion at the high binding energy side 
and its extrapolation to the low energy side (dashed line). 
ARPES data is shown in the background for the comparison 
purpose. 



self energy, it is important to have the self energy in the 
full k-o; space as is the case for our method. 

The potential importance of the above observation 
can not be understated. The largest shortcoming of 
dynamical mean field theory (DMFT) is that momen- 
tum independence of the self energy must be assumed 
(some momentum dependence can be added back in 
with cluster calculations, but this is computationally very 
expensive) 26 ! 27 . While experimental limitations forced us 
to make a few assumptions such as the background sub- 
traction, our results are encouraging that the momentum 
independent self energy assumption made in DMFT cal- 
culations of the cuprates are valid. To our knowledge, 
this is the first experimental attempt to prove the valid- 
ity of this assumption. 

Looking at ImE ER plotted in panel (b), one may won- 
der about the validity of the result near the gap region as 
ImE^ is zero while it is known that ImE Fn ^ diverges 
near the gap region 21 . The divergence of ImE Fn ^ near 
the gap is due to the fact that multiple poles of G Ful1 
(or two peak structure in A Ful1 (k,o;)) are attempted to 
be described by a single pole with a self energy. How- 
ever, by using only the occupied part, A ER effectively 
becomes a single peak spectral function. As ImE Fn ^ does 
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not diverge for a single peak spectral function, ImE ER 
calculated from an effective single peak does not neces- 
sarily diverge even though CCOC is an insulator. This 
issue will be discussed again with simulations presented 
in Fig. 3. 

In the conventional method of fitting MDCs and con- 
verting the MDC widths to EDC widths 1 , one estimates 
the width, thus ImE^, at the peak positions. We 
present the result in a similar fashion by plotting ImE^ 
along the dispersion. This provides the most accurate 
estimate of T i ER of the feature because of the possible 
momentum dependence discussed above. In Fig. 2(c), we 
plot -ImE^ along the MDC dispersion shown in panel 
(b). At this point, it is worth checking the validity of 
the new result. One way of doing so is to compare the 
ImE ER obtained by the new method with the value from 
a conventional method (if such a method can be applied). 
In the case of CCOC, a conventional way of measuring 
the half width of a quasiparticle-like peak may be appli- 
cable only to the features with the lowest binding energy. 
In the inset, we plot the EDC of the feature that has the 
lowest binding energy of 0.4 eV, hence is sharpest. The 
half width of the feature is about 150 meV. ImE ER at 
0.4 eV indeed shows a similar value. This fact gives us 
confidence in the new method. 

Looking at ImE^ curve in panel (c), one sees 
"bends" at about 0.3 and 1 eV binding energies. 0.3 eV 
is where the spectral weight of the lowest binding energy 
feature starts as can be seen from Fig. 1(a). There- 
fore, the bend at 0.3 eV may not mean any coupling 
energy. Meanwhile the bend at 1 eV is in the mid- 
dle of the dispersion. Interestingly, 1 eV binding en- 
ergy is where the kink-like dispersion appears to have 
"bosonic mode" coupling. When the kink-like feature 
from CCOC was first reported 9 , the most obvious prob- 
lem with the interpretation of it in terms of bosonic 
mode coupling was how an insulator can have bosonic 
mode coupling. Strictly speaking, an insulator can have 
bosonic mode coupling and should display kink-like fea- 
tures even though the effect would be small 28 . Even 
though it does not show the origin, our analysis firmly 
shows that ImE ER behaves as if there is a bosonic mode 
coupling at 1 eV (or 0.6 eV from the top of the band). 
This value is comparable to the energy found in doped 
HTSCsl2*H. 

In Fig. 2(d), we plot the MDC dispersion of the high 
energy features and its extrapolation as the dashed line. 
This has a strong resemblance to the LDA band and may 
serve as the bare band 9 . Also plotted in the panel is the 
£k obtained by calculating the center of the gravity of 
the spectral function. In comparing the two, one sees 
very little resemblance between them. The deviation of 
Sk from the bare band near k = 0.5A -1 (where ^ 1) is 
from the reason discussed above. Meanwhile, we believe 
the difference near k = 0, where we expect rik to be closer 
to 1, is from the fact that we cut out the spectral weight 
at high binding energies under the main valence band 29 
during the background subtraction. 
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FIG. 3: (a) A Lorentzian peak with an width of 0.3 eV rep- 
resenting a spectral function A Full (k, u). The filled area be 
low the Fermi energy represents the electron removal spec- 
tral function A ER (k,uj). (b) -ImE~ 
(dashed) of the spectral functions in ( 
peaks of 0.3 eV width, simulating a situation where there is 
a gap at the Fermi energy. The filled area again represents 
A ER (k, u). on the red solid line (blue solid line), (d) -Im£ Fu " 
(solid) and -lmT> ER (dashed) of the spectral functions in (c). 



An important question is whether the self energy de- 
fined only for the electron removal process is meaningful. 
We argue that it still bares useful information. First of 
all, the self energy defined in this way is identical to the 
usual Yi FuU when n& = 1. This is important because 
for most of the cases one analyzes the self energy for the 
occupied parts of the momentum space where nk ~ 1. 
In addition, we find that ImE^ measures the "width" of 
spectral function fairly well in various cases as illustrated 
below. Lastly, on the practical side, the method can 
also be applied when the dispersion is flat and therefore 
MDCs can not be obtained. This is very useful when the 
dispersion does not cross the Fermi energy as in insu- 
lators. To discuss if ImE ER indeed measures the width 
of the occupied part of a spectral function well, we use 
numerical simulations as illustrated in Fig. 3. The first 
case, shown in panel (a), is when the spectral weight is 
truncated by the Fermi function. From the results plot- 
ted in panel (b), we find that ImE^ still measures the 
width of the peak fairly well below the Fermi energy. We 
may also discuss cases where there is a gap at the Fermi 
energy and the spectral function splits into two due to, for 
example, charge density wave or Mott transition. To sim- 
ulate such cases, we use two Lorentzian peaks as shown 
in panel (c). From the self energies plotted in panel (d), 
we find that ImE ER is very close to the Lorentzian peak 
width while ImE Fn ^ has a large value near u=0. These 
numerical simulation results support the assertion that 
l m j]ER measures the width of the occupied part of a 
spectral feature. The case illustrated in panels (c) and 
(d) is also relevant to the discussion on the behavior of 
the self energies in the gap region. While ImE Fn ^ has a 
large value at cj=0, ImE^ stays at a finite value close 
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to the peak width. 

Our method, as it is, may not be used for multi-band 
systems. This is due to the fact that states from different 
bands may have different matrix element M(k). How- 
ever, our scheme in principle can be extended to multi- 
band systems even though it is more difficult. If it can be 
extended to multi-band system, our method will find far 
more applicability and usefulness in the analysis of self 
energies of many-body systems. This is at the present 
left for the future work. 



IV. CONCLUSION 

In conclusion, we propose using the electron removal 
self energy which is found to represent the spectral width 
very well. The method has the advantage that one can 
easily extract the self energy from a spectral function 



with arbitrary shape, allowing us to avoid making inap- 
propriate assumptions. Application of the new method 
to experimental data from CCOC reveals momentum in- 
dependence of E(k,u;) which has been an assumption in 
the conventional analysis. We expect this new method 
to be useful in strongly correlated systems where spectral 
features are often quite broad. 

Acknowledgments 

Authors would like to thank B. J. Kim for a key sug- 
gestions in programming and W. Meevasana, Jun Won 
Rhim and Jung Hoon Han for helpful discussions. This 
work is supported by the KOSEF through CSCMR and 
by the KICOS through a grant provided by MOST in 
M60602000008-06E0200-00800. Work at Los Alamos was 
performed under the auspices of US DOE. 



Electronic address: cykim@phya.yonsei.ac.kr 

1 B. A. McDougall, T. Balasubramanian, and E. Jensen, 
Phys. Rev. B 51, 13891 (1995). 

2 M. Hengsberger, D. Purdie, P. Segovia, M. Gamier, and 
Y. Baer, Phys. Rev. Lett. 83, 592 (1999). 

3 T. Valla, A. V. Fedrov, P. D. Johnson, and S. L. Hulbert, 
Phys. Rev. Lett. 83, 2085 (1999) 

4 P.V. Bogdanov, A. Lanzara, S.A. Kellar, X.J. Zhou, E.D. 
Lu, W.J. Zheng, G. Gu, J.-I. Shimoyama, K. Kishio, H. 
Ikeda, R. Yoshizaki, Z. Hussain, and Z.-X. Shen, Phys. 
Rev. Lett. 85, 2581 (2000). 

5 A. Lanzara, P. V. Bogdanov, X. J. Zhoud, S. A. Kellar, D. 
L. Feng, E. D. Lu, T. Yoshida, H. Eisaki, A. Fujimori, K. 
Kishio, J.-I. Shinoyama, T. Noda, S. Uchida, Z. Hussain, 
and Z.-X. Shen, Nature (London) 412, 510 (2001). 

6 A. Kaminski, M. Randeria, J. C. Campuzano, M. R. Nor- 
man, H. Fretwell, J. Mesot, T. Sato, T. Takahashi, and K. 
Kadowaki Phys. Rev. Lett. 86, 1070 (2001). 

7 A. D. Gromk, A. V. Fedorov, Y.-D. Chuang, J. D. Koralek, 
Y. Aiura, Y. Yamaguchi, K. Oka, Yoichi Ando, and D. S. 
Dessau, Phys. Rev. B 68, 174520 (2003). 

8 W. Meevasana, N. J. C. Ingle, D. H. Lu, J. R. Shi, F. 
Baumberger, K. M. Shen, W. S. Lee, T. Cuk, H. Eisaki, 
T. P. Devereaux, N. Nagaosa, J. Zaanen, and Z.-X. Shen, 
Phys. Rev. Lett. 96, 157003 (2006). 

9 F. Ronning, K. M. Shen, N. P. Armitage, A. Damascelli, 

D. H. Lu, Z.-X. Shen, L. L. Miller, and C. Kim, Phys. Rev. 
B 71, 094518 (2005). 

10 B. P. Xie, K. Yang, D. W. Shen, J. F. Zhao, H. W. Ou, 
J. Wei, S. Y. Gu, M. Arita, S. Qiao, H. Namatame, M. 
Tanaguchi, N. Kaneko, H. Eisaki, Z. Q. Yang, and D. L. 
Feng, Phys. Rev. Lett. 98, 147001 (2007). 

11 J. Graf, G.-H. Gweon, K. McElroy, S. Y. Zhou, C. Jozwiak, 

E. Rotenberg, A. Bill, T. Sasagawa, H. Eisaki, S. Uchida, 
H. Takagi, D.-H. Lee, and A. Lanzara, Phys. Rev. Lett. 
98, 067004 (2007). 

12 C. Kim, F. Ronning, A. Damascelli, D. L. Feng, Z.-X. 
Shen, B. O. Wells, Y. J. Kim, R. J. Birgeneau, M. A. 
Kastner, L. L. Miller, H. Eisaki and S. Uchida, Phys. Rev. 
B 65, 174516 (2002). 



K. M. Shen, F. Ronning, W. Meevasana, D. H. Lu, N. J. C. 
Ingle, F. Baumberger, W. S. Lee, L. L. Miller, Y. Kohsaka, 
M. Azuma, M. Takano, H. Takagi, and Z.-X. Shen, Phys. 
Rev. B 75, 075115 (2005). 

B. J. Kim, H. Koh, E. Rotenberg, S.-J. Oh, H. Eisaki, 
N. Motoyama, S. Uchida, T. Tohyama, S. Maekawa, Z.-X. 
Shen and C. Kim, Nat. Phys. 2, 397 (2006). 

M. R. Norman, H. Ding, H. Fretwell, M. Randeria, J. C. 

Campuzano, Phys. Rev. B 60, 7585 (1999). 

M. Z. Hasan, E. D. Isaacs, Z.-X. Shen, L.L. Miller, K. 

Tsutsui, T. Tohyama, and S. Maekawa, Science 288, 1811 

(2000). 

J. M. Luttinger, Phys. Rev. 121, 942 (1961). 
D. C. Langreth, Phys. Rev. B 1, 471 (1970). 

D. A. Shirley, Phys. Rev. B 5, 4709 (1972). 

S. Hufner, Photoelectron spectroscopy : principles and ap- 
plication, Oxford University, New York : Springer- Verlag, 
1995. 

R. Bulla, T. A. Costi and D. Vollhardt, Phys. Rev. B 64, 
045103 (2001). 

Anil Khurana, Phys. Rev. B 40, 4316 (1989); P. W. An- 
derson, Princeton University lecture notes on storonly in- 
teracting fermions (unpublished). 

M. Eschrig and M. R. Norman, Phys. Rev. B 67, 144503 
(2003). 

A. A. Kordyuk, S. V. Borisenko, A. Koitzsch, J. Fink, M. 
Knupfer, and H. Berger, Phys. Rev. B 71, 214513 (2005). 
A. A. Kordyuk, S.V. Borisenko, V. B. Zabolotnyy, J. Geek, 
M. Knupfer, J. Fink, B. Buchner, C. T. Lin, B. Keimer, H. 
Berger, A.V. Pan, Seiki Komiya, and Yoichi Ando, Phys. 
Rev. Lett. 97, 017002 (2006). 

A. Georges, G. Kotliar, W. Krauth and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (2006). 

G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 
O. Parcollet, C. A. Marianetti, Rev. Mod. Phys. 78, 865 
(2006). 

C. S. Leem, B. J. Kim, Chul Kim, S. R. Park, T. Ohta, 

E. Rotenberg, H.-D. Kim, M. K. Kim, H. J. Choi, and C. 
Kim, unpublished. 

H. Eskes and R. Eder, Phys. Rev. B 54, R14226 (1996). 



